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Spin-dependent tunneling into an empty lateral quantum dot 
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Motivated by the recent experiments of Amasha et al. [Phys. Rev. B 78, 041306(R) (2008)], 
we investigate single electron tunneling into an empty quantum dot in presence of a magnetic field. 
We numerically calculate the tunneling rate from a laterally confined, few-channel external lead 
into the lowest orbital state of a spin-orbit coupled quantum dot. We find two mechanisms leading 
to a spin-dependent tunneling rate. The first originates from different electronic g-factors in the 
lead and in the dot, and favors the tunneling into the spin ground (excited) state when the g- 
factor magnitude is larger (smaller) in the lead. The second is triggered by spin-orbit interactions 
via the opening of off-diagonal spin-tunneling channels. It systematically favors the spin excited 
state. For physical parameters corresponding to lateral GaAs/AlGaAs heterostructures and the 
experimentally reported tunneling rates, both mechanisms lead to a discrepancy of ~10% in the spin 
up vs spin down tunneling rates. We conjecture that the significantly larger discrepancy observed 
experimentally originates from the enhancement of the g-factor in laterally confined lead. 

PACS numbers: 73.40.Gk, 73.63.Kv, 72.25.Mk 
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I. INTRODUCTION 

Spintronics uses the spin, rather than the charge de- 
gree of freedom of electrons for information processing. ^'^ 
This requires the ability to create, manipulate, and detect 
spin currents and accumulations, all tasks for which semi- 
conductors seem especially promising."^ Experiments on 
few-electron quantum dots in a Zeeman field have shown 
how to manipulate two-electron spin states by elec- 
trically controlling the exchange interaction in double- 
dot systems fi— and how to convert electronic spin ori- 
entations into electrostatic voltages<^— More recently, 
Amasha et al^ reported a strong and unexpected spin 
dependence of the tunneling rate into an empty quan- 
tum dot in a Zeeman field. In these experiments a lat- 
eral quantum dot is tunnel-coupled to a narrow external 
lead. Charge sensing allows to measure the time it takes 
for an electron to enter the dot after an electric pulse has 
brought the lowest orbital level below the Fermi energy 
in the lead. Once this level is Zeeman-split the tunneling 
rates into the ground state and into the first spin-excited 
level can be extracted. Changing the magnetic field and 
the geometry of the dot, the ratio x of the tunneling 
rates for the two spin states was found to vary between 
X = 1 (symmetric tunneling) to x ^ 1 (negligible tun- 
neling into the spin excited state). It was, in particular, 
found that x ^- at large in-plane magnetic field. This 
is rather intriguing as one expects tunneling rates to de- 
pend mostly on the orbital structure of the wavefunction, 
over which a Zeeman field has no effect. Because of the 
importance that such an effect might have for quantum 
dot spintronics, these experimental results need to be 
better understood. This is one of our objectives in this 
paper. 

In III-V semiconductor heterostructures, the spin- 
orbit interactions are the usual suspect behind any 
spin-dependent effect. Several mechanisms for spin in- 
jection in tunneling structures with spin-orbit interac- 



tions have been considered, most notably based on reso- 
nant enhancement of weak spin-orbit effectSr^rJ^ spatial 
modulation of the spin-orbit interactions strength,— li^ 
spin-orbit induced mass renormalizationj^ lateral 
confinement)^^— or electron momentum filteringiiS 
Here, we extend these investigations and study electronic 
tunneling into an empty quantum dot from a single ex- 
ternal lead in presence of an external magnetic field. We 
uncover two, so far neglected mechanisms for elastic spin- 
dependent tunneling. The first one appears when tun- 
neling occurs between an extended and a confined re- 
gion with different g-factors. Tunneling being elastic, 
the splitting of the orbital states in the confined region 
determines the energy of the tunneling electron in the 
continuum. Because of the discrepancy in the g-factors, 
however, the continuum wavevectors of the tunneling 
electrons depend on the spin orientation, thus the tun- 
neling rate becomes spin-dependent. Tunneling into the 
ground-state (first spin-excited state) is favored if the 
g-factor magnitude is larger (smaller) in the lead. The 
second mechanism arises because spin-orbit interactions 
are effectively weaker in the low-energy spectrum of a 
confined region than in the continuum. The direction 
of the effective Zeeman field (the sum of the magnetic 
and spin-orbit fields) seen by the electron in the lead 
and in the dot are thus different, and this opens spin- 
non-conserving tunneling channels. Because this second 
mechanism systematically favors tunneling into the spin- 
excited state, and because our numerical investigations 
estimate a discrepancy of ~10% at most in the tunnel- 
ing rates for realistic parameters, we conclude that it 
plays only a marginal role in the experiments of Ref. |9J, 
where tunneling rates into the ground-state are system- 
atically larger and discrepancies up to almost 100% are 
observed. In contrast, we conjecture that such large spin- 
anisotropies in tunneling arise from the first mechanism 
when the g-factor in the lead is substantially larger than 
the dot and the bulk value. Strongly enhanced ^-factors 



have been observed in the lowest conduction subbands of 
quantum point contacts^"— They are usually attributed 
to the electron-electron Coulomb interactions.— We ex- 
pect that a similar enhancement occurs in narrow, later- 
ally confined leads. 

To model the experiment, we consider a two dimen- 
sional dot divided by a barrier from a semi-infinite lead of 
finite width. The potential along the lead axis is sketched 
in Figll^, with FiglTjD giving a top view of the model 
structure. We assume that tunneling is elastic, which 
has been found to be the case experimentally^^ It is easy 
to understand how a spin asymmetry in tunneling can 
arise. Consider first that there is no spin-orbit interac- 
tion. Applying an in-plane magnetic field Zeeman-splits 
electronic levels. Tunneling being elastic, the tunneling 
rate is spin independent if the g-factor is constant, be- 
cause the barrier height is the same for both spin species. 
This is sketched in Fig. [T]:. If on the other hand the Zee- 
man splitting is larger in the dot, the barrier height for 
the spin excited state becomes effectively smaller, thus 
tunneling into that state is faster than into the lowest 
spin state. The situation is depicted in Fig. [TJl. Consider 
next that spin-orbit interactions are present. They are 
effectively weaker in the dot than in the lea d^^'^^ and be- 
cause of this, tunneling can be accompanied by spin- flips 
in presence of a Zeeman field. The net result is that spin 
off-diagonal tunneling channels open up. Below we show 
that this leads to a systematically larger rate into the 
spin-excited state in the dot. This mechanism is sketched 
in Fig. [1^. 

The article is organized as follows. In Sec. II we present 
our model for the dot and lead states, and for the tun- 
neling rate. We analyze separately the t^-factor inho- 
mogeneity and the spin-orbit interactions mechanisms in 
Sees. Ill and IV, respectively. We conclude in Sec. V. In 
the Appendix we compare the tunneling formula we use 
with three common alternatives in a simplified model, to 
demonstrate its generality. 



II. MODEL 

We compute tunneling rates using the method of 
Ref. [23, which to leading order gives 



27r 



E 



(<i>dot|^y|^lcad,)|2^(^dot 



77'lcad \ 



(1) 



It describes an elastic transition from the set of lead 
states labeled by the transverse channel c, longitudinal 
wavevector fc, and spin a', into the lowest orbital state 
of the dot with spin a. The effective transition potential 
6V, is defined as the difference of the potential of the iso- 
lated dot and a dot with a nearby lead, see Fig. [T^. The 
tunneling formula, Eq. ((TJ , is discussed in the Appendix, 
where we additionally show that in one dimension it is 
equivalent to other, standardly used tunneling formulas. 
The lead and the dot wavefunctions in Eq. ([TJ are those 
of the two subsystems considered separately. We take 
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FIG. 1: (Color online) (a) Potential profile of the dot and the 
semi-infinite lead when considered separately, (b) Top view 
of the dot-lead structure, (c)-(e) Schematics of the tunneling 
mechanisms, (c) Symmetric tunneling, (d) Asymmetric tun- 
neling: a larger (smaller) g-factor in the lead enhances the 
tunneling into the ground (excited) state, (e) Asymmetric 
tunneling: spin-orbit interactions open the off-diagonal spin 
tunneling channels, which always favors tunneling into the 
excited state. 



them as eigenstates of the following Hamiltonian: 



H = k" 
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with appropriate boundary conditions. These are that 
the dot wavefunctions are zero at infinity, while the ex- 
tended states of the lead are fixed by the energy, the 
transverse channel, and the spin of the incoming wave. 
In Eq. ©, TO is the electron effective mass, k — — iV, 
B is an in-plane magnetic field forming angle ^b with 
the crystallographic x axis. It couples to electronic 
spins via the Lande 5-factor and the Bohr magneton 
/is, and cr = {(JxiCfy^az) is a vector of Pauli matrices. 
The spin-orbit interactions Hamiltonian Hso includes the 
Bychkov-Rashba, linear Dresselhaus, and cubic Dressel- 
haus terms^ 
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which we write as a momentum-dependent magnetic 
field Hso — (5/2)A*BBso(k) • a. Here, we adopt the two- 
dimensional approximation, where only the lowest sub- 
band of the perpendicular confinement is occupied. In 
this case, the in-plane magnetic fiel d has n o orbital effect 
as long as the magnetic length, ^2h/eB is larger than 
the perpendicular extension of the two-dimensional elec- 
tron gas (2DEG). This conditions is fulfilled here, since 



even at the largest magnetic fields of 7.5 Tesla we con- 
sider, the magnetic length is ~ 13 nm. Moreover, the 
strengths of the two Dresselhaus interactions are related 
by h^/2mld = 7c(fcz), with (fc^) being the quantum me- 
chanical variance of kz in the lowest subband in the per- 
pendicular confinement of the two-dimensional electron 
gas. 

Electronic confinement is described by the potential 
V{r). The semi- infinite lead is laterally confined by in- 
finite potential barriers a distance wicad apart. Close to 
the dot, it is terminated by a potential barrier linearly 
reaching a magnitude Vq over a distance L (we call this 
the "ramp"). 



l^'^"''(2:iead) = 



Vo, Xlead < 0, 

(1 - Xlcad/L)Vo, < CClead < L, 

0, L < a;icad- 



(4) 



The lead coordinate system is offset by Toffsot with re- 
spect to the dot center and rotated by an angle 0ioad 
with respect to the crystallographic coordinates x, and 
y. This is shown in Fig. [TJd. 

For the dot confinement we take an anisotropic linear 
harmonic oscillator, bounded from above by Vq in the 
lead area [see Figs.[T^ and[T}3] 
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The main axes of the dot are rotated with respect to the 
crystallographic axes by an angle 0dot- The two indepen- 
dent confinement lengths Ix and ly allow us to consider 
dot deformations as in the experiment of Ref. |9|. The 
potential Vg is used to offset the dot with respect to the 
lead. The model has no particular spatial symmetry. 

We use parameters typical for GaAs heterostructures. 
When not varied, their values are m = 0.067, 7c = 27.5 
eVA^, g = -0.39, Id = 0.63^m, kr = 2.4^m^ and a 
Fermi energy in the lead Ep — 8 meV. The geometry is 
set to qualitatively correspond to the description of the 
experimental setup in Refs. |^ and [2^, Ix = ly = 24 nm, 
corresponding to an orbital excitation energy Eoih — 2 
meV, Vo = 12 meV, L = 240 nm, (/)ioad = -7i'/4, 
'/'dot — ''■/4, and 0b — 37r/4. We offset the dot by Vg « 6 
meV to align the lead Fermi energy with the dot spin 
excited state. We offset the lead with respect to the dot 
by ^offset = 24 nm and we take wicad = 72 nm, giving two 
open transverse channels. Finally, we use a;offsGt as a free 
parameter, which we always adjust to keep the total tun- 
neling rate into the spin-excited state at 200 Hz. Below 
we vary all the parameters individually, which allows us 
to identify which are relevant for spin-dependent tunnel- 
ing. Because we consider tunneling into an empty dot, we 
neglect electron-electron interactions in the dot (beyond 
charging effects which prohibit double occupancy of the 
dot), but incorporate them in a renormalized gr- factor in 
the lead. 
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FIG. 2: Spin tunneling asymmetry due to the (/-factor inho- 
mogeneity varying (a) the (/-factor difference, (b) the length 
L of the barrier ramp, (c) the barrier height Vb, and (d) the 
dot asymmetry {l^ — ly)- When not varied, the parameters 
are ffiead = -0.39, gdot = -0.32, Vo =12 meV, L = 240 nm, 
and Ix — ly = 0. In all panels, the tunneling rate into the 
ground state has been kept fixed at 200 Hz. 



We define the asymmetry x in spin tunneling as the 
ratio of the tunneling rates into the two lowest Zeeman 
split dot states 



X = Ti/rt 



(7) 



Symmetric tunneling corresponds to x = 1, while x 
smaller (larger) than one means that tunneling into the 
ground(excited) spin state is faster. 



III. SPATIALLY DEPENDENT g-FACTOR 

As a first source of tunneling asymmetry, we consider 
the inhomogeneity in the (/-factor and neglect spin-orbit 
interactions for the time being. A slight difference of the 
(/-factor in the dot and the lead is not unexpected, as in 
semiconductor quantum dots the (/-factor depends on the 
magnetic field, state energy, or wavefunction penetration 
into the barrier along the growth direction. '^'^"■^'^ Refer- 
ence [23 reported mesoscopic fluctuations of the (/-factor 
in a quantum dot of a relative magnitude '^ 20%, see 
Fig. 3a there. 

We first discuss qualitatively how a (/-factor difference 
results in X 7^ 1. According to Eq. ([T]), the tunneling rate 
from a lead state a' to the dot state a scales with 



(8) 



T,,^, ^ i($r I*.' )r - exp(-2fc,,^,L)i(rir 

The orbital overlap depends on the wavevector 

ka'^a = ^2m[{Vo + e'-^<^) - {Eo + e^°^)]/h^, (9) 

where we split the energy levels in the dot into the orbital 
and Zeeman contributions, Eq + ejj°', and ej^/^*^ is the 



Zeeman energy in the lead. In the absence of spin-orbit 
interactions, (^'^|^'^ ) — ^aa'- If the (/-factor is the same 
everywhere, there is no tunnehng asymmetry 
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If, however, the Zeeman energies in 



the dot and the barrier differ, we get 



X ~ exp -2L 



h\V^ - E) 



(5lcad -5dot)A*sS • (11) 



A g-factor larger in magnitude in the lead gives a faster 
tunneling into the ground state, as observed in Ref. |^. 
However, assuming a dot ^-factor reduction of 20%, we 
estimate x ^ 0-95 at 7.5 Tesla, an asymmetry much 
smaller than observed. 

We numerically evaluate T „ from Eq. ^ and plot x 
in Fig. [21 Figure [5^ confirms that the tunneling into the 
ground (excited) state is preferred, if the g-factor magni- 
tude is larger (smaller) in the lead than in the dot. We 
next show in Fig.[5}D that the asymmetry is larger for tun- 
nel barriers with larger ramps. For a dot-lead ij-factor 
difference of 20%, x ~ 0-9 for ^ barrier with L k, 200 
nm, which qualitatively corresponds to the experimen- 
tal setup of Ref. [3. There is a slight dependence of x 
on Vo shown in Fig. [5J;. The trend somehow contradicts 
Eq. (jlip . because for low barriers, the dot wavefunction 
leaks out into the lead channel, and this is not captured 
by the approximations made in Eq. ([5]) . Finally, Fig. [2Ji 
shows that deforming the dot leaves x practically un- 
changed. This is expected, because the symmetry of the 
ground state depends only weakly on details of the con- 
finement potential. Most importantly, these numerical 
results on a two-dimensional model confirm the qualita- 
tive picture based on a one-dimensional toy model given 
above. For variations in the dot's (^-factor up to 20%, the 
relative tunneling asymmetry does not exceed '^ 10% and 
the symmetry of the confinement potential of the dot 
plays no role. We conclude that the reported fluctua- 
tions in the dot's (/-factor variations can not explain the 
reported large spin-asymmetry in tunneling. 

A small X requires a large argument in the exponential 
of Eq. (|lip . which can be reached for a large difference 
in Zeeman energy or kh 3> 1. The barrier width L is 
limited by the tunneling rate magnitude (200 Hz), and 
therefore cannot exceed few hundreds of nanometers. We 
just argued that variations in the dot's g-factor cannot 
account for the observed Xi if 5icad = —0.39, which is 
the standardly accepted value for the two dimensional 
electron gas in GaAs/AlGaAs heterostructures. We now 
argue that only an enhancement of the lead g-factor can 
generate a large tunneling asymmetry. A strong enhance- 
ment of the g-factor, up to \g\ =3 at magnetic field of 
16 Tesla, has been reported in Refs. [20| and [2lJ in the 
lowest conduction subbands of a quantum point contact. 
It is reasonable to expect that a similar effect occurs in 
the narrow leads of Ref. [91. Two data sets are shown 
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FIG. 3: Spin tunneling asymmetry due to the enhancement of 
the g-factor in the lead as a function of (a) the barrier length 
for g- factor values similar to those reported in Refs. [2y and 
I21I and (b) the lead g-factor for two fixed barrier lengths. 



in F ig. B^ for values of the lead 5-factor as reported in 
Ref. I20I The tunneling asymmetry is now strongly en- 
hanced and comparable to the experiment at reasonable 
barrier lengths. Alternatively, Fig. [3}3 shows the tunnel- 
ing asymmetry at two fixed barrier lengths. The (/-factor 
enhancement required for a substantial asymmetry is still 
reasonable, given that enhancement factors of up to 50 
were reported in Ref. [23. Using the enhanced g-factor is 
consistent with the fact that the tunneling in our model 
is dominated by the lowest transverse channel. 

Reference 9 showed a strong variation in x ^s the po- 
tentials defining the confinement of the dot were changed. 
It was concluded that x strongly depends on the dot's 
shape. Our theory does not account for such a depen- 
dence — the data in Fig. [2}1 show only a marginal de- 
pendence of X on I3, — Zy. There are two mechanisms by 
which changing the gate voltages defining the dot con- 
finement may influence the tunneling asymmetry. First, 
changes in those voltages can be accompanied by shifts in 
the position of the dot and thus modulations of the bar- 
rier length L and of the tunneling asymmetry."^ Second, 
shape voltages can directly influence the transverse con- 
finement in the lead tip, on which the (/-factor strongly 
depends.— 



IV. EIGENSPINOR ORIENTATION MISMATCH 

We next investigate the influence that spin-orbit in- 
teractions have on x^ assuming a homogeneous g-factor. 
The numerical results allow us to rule out the spin-orbit 
interactions as the source for the observed asymmetry, 
as the effect is too small. For the sake of completeness, 
we nevertheless explain how the spin-dependent tunnel- 
ing arises, since it may become important in materials 
with a stronger spin-orbit coupling. 



It is important to note that the effective strength of 
spin-orbit interactions is not the same in the dot and in 
the lead. In 2DEG GaAs, the spin-orbit induced mag- 
netic field is ^10 Tesla. As this is at least comparable to 
the external magnetic field, spin-orbit interactions signifi- 
cantly influence the eigenspinor orientation. On the other 
hand, spin-orbit interactions are strongly suppressed in 
the dot. The spin is almost perfectly aligned with the ex- 
ternal magnetic field, defiecting from it by a small angle 
of the order of lx,y/lso <S li ^^'^^ Due to the misaligned 
eigenspinors, spin-dependent tunneling is expected. 

We therefore consider the eigenfunctions of the Hamil- 
tonian in Eq. ([2]) for a constant potential V{r) — V. The 
ansatz^ ^(r) = exp(ik • r)^ leads to an algebraic equa- 
tion for the wavevector k and a spatially independent two 
component eigenspinor ^, 
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[B-t-B,o(k)]-rTk = 0. (12 



For given energy and wavevector, Eq. ([T2|) has two so- 
lutions, which we label by the spin index a — ±1. The 
direction of the eigenspinor ^o- is parametrized by its in- 
clination and azimuthal angles {9a,4'cr)- We call a state 
evanescent ii V > E and extended ii V < E. These 
two have purely imaginary and real wavevector k, re- 
spectively, if one neglects the Zeeman energy e^ contri- 
bution to the eigenstate total energy E. The spin-orbit 
induced magnetic field Bgo, defined by Eq. ([3]), is in- 
plane, and generates an anisotropy if the external field 
direction is varied with respect to the crystallographic 
axes. In Fig. |4^, we plot the azimuthal angle of B and 
Bgo and mark the directions where the external field is 
(anti)parallel/perpendicular with the spin-orbit field by 
the downward/upward arrows. 

For an extended spinor, Bso(k) is real. The sum of 
the spin-orbit and external fields sets the spin direction 
of the extended solutions of Eq. p^ to 



Tr/2, (f)a = 0: 



B+B^ 



+ (1 - cr)7r/2, £„ 



The Zeeman splitting energy 



(5/2)^B\/(B + B, 



(13) 



(14) 



is (minimal) maximal, if the two fields are (anti)parallel. 
This is shown in Fig.|4j:. The opposite Zeeman energy of 
the two spinors is compensated by a small change in the 
wavevectors. 

The picture changes appreciably for an evanescent 
spinor. Despite their importance for tunneling, analyt- 
ical solutions of Eq. P^ are known only for zero mag- 
netic field, ^^ while some numerical results exist for a fi- 
nite magnetic fieldj^"— A simple perturbative solution 
of Eq. ([T2|) follows if the spin-orbit magnetic field is taken 
as 



Bso(k)«B,o(ko), 



(15) 



with the unperturbed wavevector defined by fi,ko/2r7i — 
E — V < 0. Expanding the square root in Eq. (ITil) . 







a 






X 


i/l 


71- 


1 




^y 




Ml 

C 


B, 




y / 

R 










/f 


so 


^ 




0- 


^^^ r 


1 




— 1 — 




Jt/2 


20 

15 

10 

5 



b| 







V 




^^ 



d 


1 le/^ll 


-^\/ 


\/ 


- B^ t 

so 

1 1 1 


, 



7t 



2jt 



FIG. 4: Description of extended and evanescent spinors. 
(a) In-plane direction of the external and spin-orbit mag- 
netic field. At a downward/upward arrow the fields are 
(anti)parallel/perpendicular. (b) Azimuthal {<j)) and inclina- 
tion (9) angles for the spin direction of an evanescent spinor. 
[(c)-(d)] External and spin-orbit magnetic fields, and Zeeman 
energy rescaled by /x = {g/2)iJ.B for (c) extended and (d) 
evanescent spinor. 



one finds that the error of the approximation is of order 
iEso/\E - V\y - lO-'^. The spin-orbit magnetic field is 
now purely imaginary. Still, the Zeeman interaction term 
with a complex "magnetic field"— has two eigenstates 
with opposite complex eigenenergies given by Eq. (J14l) . 
whose magnitude is plotted in Fig. HJl. In the coordinate 
system with the z axis perpendicular to the plane defined 
by B and Bgo (the crystallographic coordinate system 
here), the spin direction of the evanescent eigenspinors 
is. 



+ (1 - r7)7r/2, e. 



CTC. 



(16) 



The angles 6e , and 0e are plotted in Fig. |4Jd as a func- 
tion of the external magnetic field orientation. The two 
spinors are orthogonal and in-plane only if the two mag- 
netic fields are (anti)parallel (downward arrows). The 
out-of-plane spinor component is maximal if the two 
fields are orthogonal (upward arrows). 

With the Zeeman energy given by Eqs. ([Ti)) and (|15l) . 
the wavevector k is obtained from the quadratic Eq. (|12l) . 
The imaginary part of the Zeeman term and the kinetic 
energy cancel such that the total energy E is real. To 
illustrate the solutions, we consider the wavevector fixed 
along ko so we can write it as k = ( 1 — Kz / fco — T-Kso / ^o ) ko 
introducing two real parameters Kz/so- To leading order, 
these are given by 



= ± 



2me 



(17) 



in particular, they are opposite for the two spinors. 

Consider now that a coherent superposition of the two 
eigenspinors propagates inside the barrier. The real part 
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FIG. 5: Influence of the magnetic field strength and orien- 
tation on the spin-orbit induced spin tunneling asymmetry, 
(a) At a downward/upward arrow the external and spin-orbit 
magnetic fields are (anti)parallel/perpendicular. (b) The 
asymmetry if only linear (lin), cubic (D3), or both (dashed) 
spin-orbit interactions are present, (c) Asymmetry for varied 
out-of-plane field orientation. This is the only case in this 
manuscript, where the vector potential A contributes to the 
kinematic momentum operator, Kk = — iftV -f eA. (d) Asym- 
metry as a function of the in-plane field strength. 



of the Zeeman energy leads to different decay lengths for 
the two components, while the imaginary part leads to a 
coherent rotation of the initial spin orientation. If both 
the magnetic field and the spin-orbit interactions are 
present, the spin will make a spiral, both decaying and 
rotating at the same time, over the length scales 1/kz, 
and 1/kso, respectively. The evanescent states show an 
additional possibility for the spin selection through spin 
dependent penetration lengths. Note, however, that both 
eigenspinors are additionally damped over the length 
1/fco, which is usually much shorter than 1/kz,so, strongly 
limiting the achievable differentiation for the two spin 
species. 

Before we present numerical results we note an inter- 
esting singular limit when the external and the spin-orbit 
magnetic fields are orthogonal and of the same magni- 
tude. In this case the out-of-plane spinor misalignment 
is maximal and both spinors are oriented along the same 
direction (the z axis). 

We analyze our numerical results for the tunneling 
asymmetry using the just developed picture of the ex- 
tended and evanescent spinors. In Fig. [S^ we vary the 
direction of the external magnetic field. If the two mag- 
netic fields are (anti)parallel (downward arrows), the spin 
quantization axis throughout the structure is the same 
and only spin-diagonal tunneling is possible. The asym- 
metry arises from the spin-orbit contribution to the Zee- 
man energy in the lead, which gives rise to spin depen- 
dent discrepancies in the wavevector. When the mag- 
netic fields are not parallel, the eigenspinor orientations 
in the lead, under the barrier and in the dot are all dif- 
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FIG. 6: Dependence of the spin-orbit induced spin tunneling 
asymmetry on (a) the barrier height Vb, (b) the barrier ramp 
L [see Fig. [T^], (c) the dot-lead offset yoffsct [see Figs.[T^ and 
[T]d], (d) the lead width TO„idth. The widths at which the lowest 
three transverse channels become open are indicated. In all 
panels, the tunneling rate into the ground state has been kept 
fixed at 200 Hz. 



ferent, which opens off-diagonal spin tunneling channels. 
The spin dependent tunneling arises from different effec- 
tive potential barriers in these off-diagonal channels. The 
wavevectors in Eq. ^ fulfill [note that a ^ a' refers to 
an eigenspinor direction in the lead {a) and the dot {(j')^ 
which are not necessarily along the same axis] 



/c4,^t > ^ 



't^t 



H^i 



>fc. 



t^;- 



(18) 



They are shown as arrows of different lengths in Fig. [T^. 
The tunneling into the excited state is always preferred 
and the effect is maximal if the spin-orbit and the mag- 
netic fields in the lead are perpendicular. 

Figure [SJs shows that the asymmetry is dominated by 
the linear spin-orbit terms, with the cubic Dresselhaus 
contribution smaller by a factor of ~5. We checked (but 
do not show) that the contributions from different spin- 
orbit terms are additive with good accuracy. In Fig. [5j: 
we consider an out-of-plane magnetic field and show that, 
despite the presence of orbital effects, the asymmetry is 
of similar magnitude as for an in-plane field. Finally, we 
show in Fig. [SJi that the asymmetry grows with the mag- 
netic field. This is expected, as a larger Zeeman energy 
makes the difference of wavevectors in Eq. (TTSl) larger. 

Figure [H] shows the dependence of the effect on the 
barrier potential. The tunneling asymmetry is more pro- 
nounced for lower (panel a) and longer (panel b) barriers 
because the longer the barrier, the stronger the suppres- 
sion of the spin down-to-up tunneling channel relative 
to the up-to-down channel. We next see in Fig|BJ: that 
the dot's symmetry has little importance since varying 
the dot-lead offset changes the tunneling asymmetry only 
marginally. Similarly, Fig. |BJi shows that the lead width 
is important only around the pinch-off point, where the 



upper spin state disappears from the lead. Because the 
tunnehng is dominated by the lowest transverse channel 
contribution, opening additional channels has a hardly 
noticeable effect on x- 

Our numerical data therefore show that spin-orbit in- 
duced spin tunneling asymmetries reflect evanescent and 
extended spinor solutions of Eq. ([T2|) but that the effect is 
too small to account for experimental observations. Here 
also, the dot symmetry plays only a minor role. 



The fields of up to 7-8 Tesla used in Ref. ^ however corre- 
spond to a magnetic length significantly larger than the 
typical thickness of the 2DEG . Possible weak diamag- 
netic effects, such as a slight shift or a compression of 
the 2DEG, seem to affect both the lead and the dot in 
the same way and should not result in the spin tunneling 
asymmetry. 
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V. CONCLUSIONS 

We investigated possible spin dependencies in the tun- 
neling of an electron into an empty lateral quantum dot 
in presence of a magnetic field. We used a general two di- 
mensional model and examined the parameters' space in 
great details, by varying the tunnel barrier geometry, the 
spin-orbit interactions strengths, the magnetic field ori- 
entation and strength, the dot confinement potential, the 
dot 5-factor, and the dot vs lead orientation. We found 
that the spin tunneling asymmetry does not exceed 10% 
for realistic parameters, much too small compared to the 
results of Ref. 9'. We therefore postulate that the ob- 
served asymmetry originates in an anomalously enhanced 
lead ^-factor, because of the lateral confinement in the 
lead. Such an enhancement has been reported in Ref s. 1201 
andil 

We expect the error stemming from approximations in 
our model to be negligible compared to the specificities 
of the tunnel barrier profile, about which not much is 
known in gated quantum dots. On the other hand, we 
do not expect the asymmetry to arise because of specific 
barrier shape, as seemingly related effects are present in 
different structures4S 

Our findings provide cues for further experimental 
checks of the origin of the spin tunneling asymmetry. As- 
suming the g-factor mechanism is at play, the tunneling 
asymmetry growth is accompanied by a reduction in the 
tunneling rates themselfs. On the other hand, the spin- 
orbit interactions induce anisotropy with respect to the 
crystallographic directions and we identify specific direc- 
tions of the magnetic field, which give maximal/minimal 
spin tunneling asymmetry. This mechanism might be of 
relevance in other materials, such as InAs4i We note, 
but do not show, that our numerics indicate that these 
two effects are essentially additive (independent). 

Among alternative explanations for the observed asym- 
metry, one could consider that single electron tunneling 
is dynamically assisted by phonons and/or nuclear spins. 
However, one does not expect these to induce a spin- 
dependence, as phonons do not couple and nuclear spins 
couple only very weakly to the electronic spin. The de- 
pendence may arise due to an inhomogeneous collective 
nuclear field, which would be similar to, and presumably 
indistinguishable from, the inhomogeneous g-factor that 
we studied here. We note finally, that the finite extension 
of the 2DEG in the z-direction allows for orbital effects. 
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Appendix A: Comparison of various formulas for 
tunneling 

We compare four different formulas for the dot to lead 
tunneling rate in one dimension. We differentiate scat- 
tering approaches, based on eigenstates of the full sys- 
tem Hamiltonian from perturbative approaches, where 
the dot and the lead are considered separately. We find 
that the perturbative approach used in the main text 
agrees very well with scattering approaches. A similar 
correspondence is expected for tunneling into a quasi 
one-dimensional channels of a finite- width lead which we 
consider in the main text. 

We set up a one dimensional square barrier model with 
the Hamiltonian H = h^k^/2m + V{x). The potential 
V{x), depicted in Fig. [7^, describes a dot of width ^o, 
separated from a semi-infinite lead by a barrier of width 
L. The dot is offset with respect to the lead by Vg and the 
barrier height is Vq. We compare the rates for tunneling 
to all quasi-bound states of the structure. 

We take the Hilbert space to be spanned by a bound 
state \(f>o) and a set of (orthonormal) extended states 
{\'4'n)}n^o- We expand the time dependent system wave- 
function as^— 



I*) = e-'"°*ao|0o) + J2 e"'""*a"|^r, 



(Al) 



introducing so far unspecified frequencies uji and time 
dependent coefficients ai{t). The time dependent 
Schrodinger equation for |^) reads 

ihdtao = ~{(j)Q\hujQ - H\(j)Q)aQ (A2a) 

n 

ihdtan = -e''^"°\^Pn\ifidt + hwo - H\cj)o)ao (A2b) 



where we denote uJij = u)i — ujj . 









FIG. 7: (a) Full scattering eigenstate (thick line) and the con- 
finement potential (thin solid line) of the dot-barrier-lead sys- 
tem, (b) Potential of an isolated dot and its bound eigenstate 
and the perturbation 5V. (c) Extended eigenstate for the 
same problem, (d) Potential profile and an extended eigen- 
state of an isolated lead. 
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FIG. 8: Comparison of two scattering approaches, (a) Wigner 
method: the refiection phase with a magnified resonance, (b) 
The flux method: lifetime according to Eq. (|A10|l . The meth- 
ods' schematics are on the right. The model parameters [de- 
picted in Fig. [Si] are: Vg = 7 meV, Vo = 13 meV, L = 50 
nm, and Iq = 100 nm. 



Let us now consider the lead as a perturbation of an 
isolated dot, as shown in Fig. [71d, that is, H = Hq + 5V , 
where 

5V = -9{x -lo- L)6V = -e{x - /o - L)Vo. (A3) 

Accordingly, we choose |0o} and \^pn) = IV-"^"') as the 
bound and extended eigenstates of Hq, depicted in 
Figs.[7}D and[7J:, respectively, and hoJi their energies. The 
formula for the tunneling rate follows from Eq. (jA2p as^ 

Ti = T E \{<Po\SV\i^^^°')fS{h^o S^r + 5V), (A4) 



if one neglects the transitions between extended states 
[m j^ n terms in Eq. (jA2bp ] and the bound state energy 

shift, i.e. ((l3o\SV\(l}o) « 0. 

Equation (jA4[) turns out to be a poor estimate. 
Gurvitz et alM-i^ suggested to improve it by replacing 
the extended states of the isolated dot |V'5?°*) [see Fig. [71;] 
by the extended states of the isolated lead j'^Jj^'^^) [see 
Fig. Ell]. We obtain Eq. (H]), 



27r 



El 



\{Msv\4>T'')rs{f^o~hw, 



lcad\ 



(A5) 



Equation (jA5l) is equivalent to the Bardeen formula for 
tunneling,'*'* which comes from the following derivation. 
In Eq. (|A1[) . we choose \(J)q) and IV'n) to be the dot bound 
state [Fig.jTt)] and lead extended states [Fig.[7li]. Transi- 
tions stem from those states being neither orthogonal nor 
eigenstates of H. Neglecting continuum to continuum 
transitions again in Eq. (jA2bl) and the time derivatives 
of the coefficients a^ on the right hand side of Eq. (|A2[) 
and using that for uq — aj„ we can write 



/)o|-ff - hWr\i}n) = ((/-ol^ - ^olV'n) = (00 It^V^j V'n) • 



(A6) 



We then get Eq. (|A5|l from Eq. ((A2|) . 

We next turn our attention to scattering approaches. 
The first one is due to Wigner.— We write a scattering 
state (an unnormalized eigenstate of iJ) in the lead as 



(a;ioad|^tuii) = e 



'\kx\. 



'\hx\r 



'\ip 



(A7) 



Current conservation requires that the reflection results 
only in a phase shift ^p. Consider now a wave packet 
of energy E and wavevector — fc (velocity v — hk/m), 
scattering at the barrier from the right. The reflected 
wave phase shift in Eq. (|A7p is interpreted as an action 
phase due to time of the tunneling ip — Et/h = 2ET^^ /h. 
Differentiating with respect to the energy we get 



dip _ 2Edr~^ 
'dE ~ IT dE 



2r- 



(A8) 



The quasi-bound state lifetime F"* is sharply peaked 
around the resonance energy, where its first derivative 
is zero. The tunneling rate follows as 



2 l'dip_ 
h KdE 



(A9) 



The reflection phase ip is plotted in Fig. [S^. The sharp 
resonances where the phase jumps by 27r are evident. 

The second scattering approach we consider is the fol- 
lowing. Take a steady state where a constant flux im- 
pinges on the barrier. A part of it is directly reflected 
while a part tunnels through. The latter (influx) in- 
creases the probability accumulated behind the barrier, 
compensated by a constant leak out (outflux). The tun- 
neling rate can be deflned as the ratio of the influx and 
the accumulated probability,-^ 



i-i = ^r — 



t'a|*Mr"'(2^1cad=0)|2 



Q L 



dot+barricr 



|*fuii(a;)|2da; 
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FIG. 9: Comparison of tunneling formulas varying the width 
of the potential welHo ■ Parameters are as in Fig. [S] Tunneling 
rate for the dot bound states computed using (a) Eq. (|A5|) 
and (b) Eq. (|A4|l . where only two states are shown for clarity, 
(c) Comparison of Eqs. (|A5|I and HA9|) . (d) Comparison of 
Eqs. (|A5|) and (|A9|) . with the latter evaluated at the energy 
of the isolated dot state. In (c) and (d) , the relative difference 
SVij = |r, - Vj\/Vi is plotted. 



The coefficient a quantifies how much of the left-going 
probabiUty flux tunnels through (a part is reflected right 
at the barrier). In analogy with an infinite barrier, we 
suppose the phase of a refiected wave changes sign at 
the reflection point, which gives a = cos^((p/2). As on 



resonance f = 7r/2, we get a = 1/2. Figure [8)d shows 
the quasi-bound state lifetime computed by Eq. (jAlOp . 
The resonant energies correspond to those obtained in 
the Wigner method. 

We compare the results of the discussed approaches in 
Fig. [9] We take the Bardeen-Gurvitz formula, Eq. (jASp . 
as reference and plot the corresponding rates in Fig. |9K- 
The Fermi Golden rule of Eq. (jA4l) shows spurious reso- 
nances and can differ by orders of magnitude, as shown 
in Fig. ^\p. As Ty, = r4 within numerical precision of our 
software, we compare only the Wigner formula to the 
reference in Fig. [SJ: and find an excellent agreement even 
improving as the resonance energy is reduced. 

To investigate the spin tunneling asymmetry, we have 
chosen to use the Bardeen-Gurvitz formula, Eq. (jA5l) . 
Out of those we discussed, it presents the advantage that 
it is based on the eigenstates of the isolated subsystems, 
which are much easier to obtain than the scattering eigen- 
states of the full system. An additional disadvantage of 
the scattering approaches is that as the tunneling rate 
drops, the resonances become narrower and harder to 
spot. Therefore, instead of searching for the them nu- 
merically, one can take the energy of a dot bound state as 
the approximation for the resonance position. In Fig. [8ji 
we evaluate the rate in the Wigner approach at such ap- 
proximate resonance energy. For states well below the 
barrier, we still find a very good correspondence with 
the Bardeen-Gurvitz formula. 
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